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In a recent paper by Link, it was pointed out that the the standard picture of the neutron 
star core, composed of a mixture of a neutron superfluid and a proton type-II superconductor, is 
inconsistent with observations of long period precession in isolated pulsars. In the following we 
will show that intervortex force between the magnetic flux tubes may be attractive resulting in a 
type-I (rather than type-II) superconductor. In this case the magnetic field cannot exist in the form 
of magnetic flux tubes, supporting Link's observation. This behavior of the system is due to the 
strong interaction between the proton-neutron Cooper pairs, which was previously ignored. We also 
calculate the critical magnetic fields H c and H C 2 for type-I/II superconductors. These results also 
support our claim of type-I superconductivity in the cores of neutron stars. 

PACS numbers: 



I. INTRODUCTION 

Recently, it was pointed out by Link pj that obser- 
vations of long period precession [2| may be inconsistent 
with the standard picture of the interior of a neutron star. 
In the conventional picture, the extremely dense interior 
is mainly composed of neutrons, with a small amount of 
protons and electrons in beta equilibrium. The neutrons 
form 3 P2 Cooper pairs and Bose condense to a super- 
fluid state, while the protons form 1 Sq Cooper pairs and 
Bose condense to give a superconductor as well (see e.g. 
Ref . [2j for a recent review) . It is generally assumed that 
the proton superconductor is a type-II superconductor, 
which means that it supports a stable lattice of mag- 
netic flux tubes in the presence of a magnetic field. This 
belief is based on simple estimations of the coherence 
length and the London penetration depth which ambigu- 
ously imply a type-II superconductivity. In addition, the 
rotation of a neutron star causes a lattice of quantized 
vortices to form in the superfluid neutron state, similar 
to the observed vortices that form when superfluid He is 
rotated fast enough. The axis of rotation and the axis 
the magnetic field are not aligned. This, coupled with 
the fact that the two different types of vortices interact 
quite strongly lead to the suggestion that the observed 
long period precession may require reexamining the pic- 
ture of type-II superconductivity inside neutron stars [lj 
that follows from the standard analysis when only a sin- 
gle proton field is considered. A possible resolution of 
this "apparent contradiction" was suggested in our recent 
letter Q , where we argued that if one takes into account 
that the Cooper pairs of neutrons are also present in the 
system and they strongly interact with proton Cooper 
pairs, the superconductor may in fact be type-I, even 
when a naive analysis that only includes the proton de- 
grees of freedom seems to indicate type-II behavior. If 
this scenario is realized in nature, this means that the 
interior of neutron stars would exhibit the Meissner ef- 
fect and therefore would not support a stable lattice of 



magnetic flux tubes. This would resolve the apparent 
discrepancy between the observation of long period 
precession 2] and the typical parameters of the neutron 
stars which naively suggest type-II superconductivity in 
neutron stars. The main goal of this paper is to provide 
a detailed analysis with complete calculations supporting 
the claim of our short letter Q. 

It is well known that type-II superconductors have 
magnetic flux tubes in the presence of a magnetic field. 
In the interior of neutron stars, which is a mixture of 
neutron and proton superfluids, the proton superfluid 
will form a vortex lattice of magnetic flux tubes if the 
superconductor is type-II. Inside the core of these vor- 
tices, the proton condensate vanishes, and the core is 
filled with normal protons resulting in the restoration of 
the broken U(1)em symmetry. For accepted estimates of 
proton correlation length and London penetration depth, 
the distant proton vortices repel each other leading to 
formation of a stable vortex lattice. This is the standard 
picture realized in conventional type-II superconductors. 
However, there are many situations where this picture 
will be qualitatively modified. For example, if there is a 
second component (such as a neutron component in our 
specific case), it may be energetically favorable for the 
cores of vortices to be filled with a nonzero condensate 
of this second component, as it was originally suggested 
in the cosmological context by Witten y|, |(| . There are 
numerous examples of physical systems where this phe- 
nomena occurs: superconducting cosmic strings in cos- 
mology, magnetic flux tubes in the high T c supercon- 
ductors, Bose-Einstein condensates, superfluid 3 He, and 
high baryon density quark matter [1 la. R. H H fToL ITU . 
Given this, one might guess that such nontrivial vortex 
structure may occur in the core of neutron stars where 
we have another example of a two component system. 
We shall argue in what follows that if the interaction be- 
tween proton and neutron Cooper pairs is approximately 
equal (a precise condition of this "approximately" will 
be derived below), the vortex- vortex interaction will be 
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modified and the system will be a type-I superconductor 
where the magnetic field is completely expelled from the 
bulk. 1 The main assumption that we are making is that 
the interactions between the proton and neutron Cooper 
pairs are approximately equal, leading to an approximate 
C/(2) symmetry. We believe that this assumption is jus- 
tified by the original isospin SU(2) symmetry of the neu- 
trons and protons. The result of this is that the proton 
vortices or magnetic flux tubes have nontrivial core struc- 
ture. The superfluid density of the neutrons is larger in 
the vortex core than at spatial infinity. In addition, the 
size of the vortex core and the asymptotic behavior of 
the proton condensate is modified due to the additional 
neutron condensate. The most important result of these 
effects is that the interaction between distant proton vor- 
tices may be attractive in a physically realizable region of 
parameter space leading to type-I behavior: destruction 
of the proton vortex lattice and expulsion of the magnetic 
flux from the superconducting region of the neutron star. 

This paper is organized as follows. In Sec. II we will 
introduce the free energy that describes a two compo- 
nent (neutron- proton) supcrfluid/supcrconducting sys- 
tem and show that the addition of the second component 
(neutron) leads to nontrivial core structure, which alters 
the properties of the magnetic flux tube. In Sec. Ill, we 
will give two different calculations of the interaction be- 
tween two widely separated vortices. These two different 
methods lead to the same conclusion, that the interior of 
a neutron star may be a type-I superconductor. In the 
subsequent section, Sec. IV, we will calculate the critical 
magnetic fields H c and H C 2 (the magnetic fields above 
which superconductivity is destroyed in type-I and type- 
II superconductors, respectively). These results confirm 
our findings in Sec. Ill that type-I superconductivity may 
occur in the interior of neutron stars. In Sec. V we will 
end with concluding remarks on possible implications of 
our results. Specifically, we will comment on the possible 
nature of glitches (observed in many systems) when the 
environment of the neutron star core is a type-I super- 
conductor. 



II. STRUCTURE OF MAGNETIC FLUX TUBES 

We start by considering the following effective Landau- 
Ginsburg free energy that describes a two component su- 
perfluid Bose condensed system. In our system, we have 
a proton condensate described by ip± and a neutron con- 
densate described by ip2- We do not consider the normal 
component of the protons and neutrons with their spe- 
cific excitations, only the superfluid component. The tpi 



In reality, the magnetic field must be present in the neutron 
star interior. The way this picture may be realized in nature is 
through the formation of domains of superconductor matter and 
normal matter. 



field with electric charge q (which is actually twice the 
fundamental charge of the proton, q — 2|e|) interacts 
with the gauge field A, with B = V x A. The two di- 
mensional free energy reads (we neglect the dependence 
on third direction along the vortex such that T measures 
the free energy per unit length): 

+ U+m+iiM*!! 2 )) (i) 

where m c = 2m and m is the mass of the nucleon. Here 
we have moved the effective mass difference of the proton 
and neutron Cooper pairs onto the interaction potential 
V. In the free energy given above, we have ignored the 
term coupling the proton and neutron superfluid veloci- 
ties, which gives rise to the Andreev-Bashkin effect |l2j . 
as it is not important in our discussion. We have also 
ignored the fact that the neutron condensate has a non- 
trivial 3 P2 order parameter as only the magnitude of the 
neutron condensate is relevant to the effect described be- 
low. 

The free energy Q is invariant under a U(l)i x U(l)2 
symmetry associated with respective phase rotations of 
fields ipi and ■02, which corresponds to the conservation 
of the number of Cooper pairs for each species of parti- 
cles. Moreover, we assume, that the proton and neutron 
Cooper pairs interact approximately in the same way. 
Therefore, the interaction potential V can be approxi- 
mately written as ^(|V>i| 2 , IV^I 2 ) « E/(|-0i| 2 + l^l 2 )- 

Naively, one might think that such an assumption 
can not be justified due to the well-known differences in 
structure of the Fermi-surfaces and gaps (which strongly 
depend on the Fermi energies) of neutron and pro- 
ton Cooper pairs in neutron stars. However, this does 
not imply that the interaction at large distances (much 
larger than the inverse gap) between proton and neutron 
Cooper pairs (not between protons and neutrons) is very 
different or that their respective Bose chemical potentials 
/ii (do not confuse with original Fermi chemical potentials 
corresponding to protons and neutrons) are very differ- 
ent. We expect that at very large distances which are 
relevant for the analysis of the vortex-vortex interaction, 
the internal structure of the gap as well as the differences 
in the densities of proton and neutron Cooper pairs n\^. 
(do not confuse with proton and neutron densities) are 
not very important provided that scattering lengths of 
different Cooper pair species are approximately the same. 
In the case where nid = m u and electromagnetic interac- 
tions are neglected, we expect that the scattering lengths 
for proton and neutron Cooper pairs would be exactly the 
same. When md ^ m u we expect that the effect of the 
asymmetry due to the renormalization with the environ- 
ment would be proportional to — m u . In principle, 
some mild singularity such as ln(rrid — m u ) may occur 
due to the renormalization. However, we do not expect 
that a strong singularity (such as (md — Triu)^ 1 ) that is 
capable of eliminating the original small factor (md — m u ) 
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will appear. In other words, if interaction of a Cooper 
pair at k — > is symmetric for protons and neutrons, it 
is not essential that the system itself is an asymmetric 
one, with n\ ^ n 2 - Such an asymmetry can be easily ad- 
justed by slightly different Bose chemical potentials, see 
below, such that the effective Lagrangian remains sym- 
metric. Therefore, when the potential is expressed in 
terms of Bose chemical potentials (macrocanonical de- 
scription) rather than in terms of densities, we expect 
the potential retains its original symmetries. 

Therefore, we believe that our approximate U(2) sym- 
metry is somewhat justified by the original isospin sym- 
metry of the original protons and neutrons, however this 
symmetry is not exactly equivalent to the conventional 
isotopical SU(2) symmetry. In reality this U(2) sym- 
metry is explicitly slightly broken, and the potential V 
has a minimum at \ipi\ 2 — ni, |^>2| — n 2, where the 
bulk proton and neutron superfluid densities ri\ and ri2 
are both non-zero. Hence in the ground state, (I^Ail 2 ) — 
7ii, i = 1,2, and both £7(1) symmetries are spontaneously 
broken. In our letter we presented some general ar- 
guments for a generic potential Vd^il 2 , 1^2 1 2 ) which has 
the symmetry properties discussed above. In order to 
be more specific and give the details of our calculations, 



in this work we will use the following standard 
potential: 
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where fi^ is the chemical potential of the i component 
and dij is proportional to the scattering length Zy be- 



tween the i and j th components, ay = AttH lij/m c . 
In our nonrelativistic formalism, the fields ipi have en- 
ergy dimension 3/2 and therefore the particle density is 
ni = {\iPi\ 2 ),i = 1,2. 

Let's analyze the vacuum structure of our potential 
(J2J. In the limit of exact 17(2) symmetry fii = fi, 
fly = a, the vacuum manifold is given by the three 
sphere \ipi\ 2 + IV^I 2 = fi/a. We are, however, interested 
in the case when the U(2) symmetry is explicitly bro- 
ken to U(l) x U(l), giving a particular pattern of proton 
and neutron condensation (IV'il 2 ) — n i- It is natural that 
even very small £7(2) symmetry violating terms are ca- 
pable of selecting a particular vacuum on the original 
degenerate manifold. First, consider the case when the 
fourth order couplings are fully degenerate ay = a, while 
the chemical potentials are slightly different fi\ = fi — 8fi, 
fi2 = fi + 8fi. In this case, the condensation pattern 
is determined solely by the sign of Sfi. If Sfi > then 
neutrons condense, (IV^I 2 ) = ^2/^22, while protons re- 
main uncondensed, (ipi) = 0; if Sfi < then protons 
condense, (IV'il 2 ) — Mi/ a iii while neutrons remain un- 
condensed, (1P2) = 0. Observe, that a very small £7(2) 
violating change in chemical potentials fit, /i 2 produces 
a very large asymmetry of proton and neutron Cooper 
pair densities m, n 2 . 



Now consider a more general situation in which all 
chemical potentials fit,fi2, and fourth order couplings 
an, ^22 7 0-12 are non-degenerate. In this case, a new phase 
is possible, where both proton and neutron condensates 
appear. In this phase, the proton and neutron Cooper 
pair densities are given by: 



»2 



Q22M1 ~ a 12^2 

2 ; 
ail022 — dl2 

aufJ-2 - Q12M1 
a? 



(3) 
(4) 



0-110-22 — "12 

This particular vacuum will be realized if and only if, 

a 22 Mi - 012^2 > 0, (5) 
«nM2 ~ ai2Mi > 0. (6) 

Notice that conditions ll- r )l(ill imply an CI22 — 0,12 > 0. If 
conditions il~>l(il) are not satisfied then only one conden- 
sate will appear, similar to the case of degenerate ay's 
already described above: if /i|/an < /i| /a 2 2 then n 2 = 
/i 2 /a 22 , ni = 0; if ^ 2 /an > ^\l a 22 then n\ = Hi/an, 
n 2 = 0. 

Throughout the rest of the paper, we will be working in 
the sector where both ipi and i\)2 obtain non-zero expecta- 
tion values, as this is the situation, which is believed to be 
realized in neutron star interiors. In this case, the chemi- 
cal potentials fit, \i2 are fixed by the equilibrium Cooper 
pair densities ri\, n 2 through equations fl.'il 11) . It is conve- 
nient to assume the following particular parametrization 
of explicit U(2) violation: /ii = fj, — 8fi, /i 2 = fi + Sfi, 
where Sfi/ fi <C 1, and an = 022 = a, a 12 = a — 8a, where 
8a/ a <^ 1. In this case the equations for Cooper pair 
densities ll.il II) reduce to: 
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(7) 
(8) 



and conditions for stability of the vacuum lloKil) reduce 
to 



\Sfi\ 8a 8a 
fi 2a — 8a 2a 



(9) 



The approximation made in the second part of equations 
(|7I8I9|I neglects terms of order 8a/ a, which are small in 
the limit of approximate J7(2) symmetry. An important 
quantity for the analysis that follows will be the ratio 
of proton Cooper pair density to neutron Cooper pair 
density, 7 = n\/ni. A typical value of 7 in the core of 
a neutron star is expected to be quite small, 5 — 15%. 
Thus, we will often use the limit 7 -C 1 in our discussion. 
However, one should remark here that our qualitative re- 
sults do not depend on the value of 7, as we will explain 
in what follows. As already mentioned, the strong devi- 
ation of 7 from 1 does not imply a large asymmetry in 
the interaction between different species of particles. As 
is clear from equations 1718191) a value of ") very different 
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from 1 can be achieved by small U{2) violating terms 
proportional to da, Sfj, in the free energy. 

The Landau-Ginzburg equations of motion following 
from the free energy are: 



in the standard fashion: 



2m c ^ Tic ^ ^ X 



+ (a - Sa)^ 2 ^, (10) 



V 2 V>2 = -M2^2 + a\t/j 2 \ 3 + (a - ^IV'il 2 ^, (11) 



2m_ 
V x (V x A 

47T 



-^[^•(V-^A)^-M (12) 
2m c c nc 



Now let's investigate the structure of proton vortices, 
which exist due to the spontaneous breaking of the £7(1) i 
symmetry. Such vortices are characterized by the phase 
of the -01 field varying by an integer multiple of 2ir as 
one traverses a contour around the core of the vortex. 
By continuity, the field must vanish in the center of 
the vortex core. Up to this point, it has been assumed 
that the neutron order parameter ip2 will remain at its 
vacuum expectation value in the vicinity of the proton 
vortex. As we have already remarked, this is not the case 
in many similar systems. Actually, it can be shown by us- 
ing the equations of motion obtained from the free energy 
Q that for most potentials V, which are approximately 
invariant under the U(2) symmetry, it is impossible for 
the ip2 field to remain constant when the ipi field varies 
in space. On the other hand, from the energetic point of 
view, given that the potential V can be written approx- 
imately &s V ~ U(\tpi\ 2 + l^hl 2 ), one can argue (in part 
based on previous work H II 0, II H E El) that it is 
favorable for the ip2 field to increase its magnitude in the 
vortex core to compensate the decrease in the magnitude 
of tpi. 

So, anticipating a non-trivial behavior of the neutron 
field -02) let's adopt the following cylindrically symmetric 
ansatz for the fields describing a proton vortex with a 
unit winding number: 



1p2 



f(r) e ie , 
V™2 g(r), 

ftca(r) - 



(13) 



where (r, 8) are the standard polar coordinates. Here we 
assume that the proton vortex is sufficiently far from any 
rotational neutron vortices, so that any variation of ip2 
is solely due to the proton vortex. The functions /, g, 
and a obey the following boundary conditions: /(0) = 0, 
/(oo) = 1, ff'(0) = 0, g(oo) = 1, o(Q) = 0, and a(oo) = 1. 
We see that the fields i\>\ and ip2 approach their vacuum 
expectation values at r = oo. 

The London penetration depth A and the coherence 
length £ of the proton superconductor will be introduced 



Airq 2 n\ ' 

/ h 2 

2m c n\a 



(14) 
(15) 



We wish to find the asymptotic behavior of fields ipi, 
ip2 and A far from the proton vortex core, as this will 
determine whether distant vortices repel or attract each 
other. The asymptotic behavior can be found analyti- 
cally by expanding the fields defined in l|13l) : 



f(r) 
9{r) 
o(r) 



1 + F(r), 
1 + G(r), 
l-r£(r), 



(16) 



so that far away from the vortex core, F, G,rS <C 1 and 
F, G, S — * as r — > oo. This allows us to linearize far 
from the vortex core the equations of motion (|1(JI11I12|) 
corresponding to the free energy 0]) to obtain: 



dr 2 
S" 



ld_ 

r dr 
1 



M 



(17) 
(18) 



where all derivatives are with respect to r and the matrix 
M mixing the fields F and G is, 



M 



4m c 



a a — 5a 
a — 5a a 



Here we assume that (rS) 2 <C F, G, i.e. the supercon- 
ductor is not in the strong type-II regime (this is justi- 
fied since we are only attempting to find the boundary 
between type-I and type-II superconductivity). The so- 
lution to Eq. I |18| l is known to be: 



Ca 
A 



Ki(r/\) 



(20) 



where K\ is the modified Bessel function and Ca is an 
arbitrary constant. The remaining equation i|17|) can be 
solved by diagonalizing the mixing matrix M. In previ- 
ous works the influence of the neutron condensate on the 
proton vortex was neglected, which formally amounts to 
setting the off-diagonal term Mi 2 in Eq. 119|) to 0. In 
that case, one can assume that the neutron field remains 
at its vacuum expectation value, i.e. G = 0, to obtain, 



F = C F K (V2r/0 



(21) 



where Kq is the modified Bessel function. It is esti- 
mated that A ~ 80 fm and £ ~ 30 fm 0, which leads to 
k = A/£ ~ 3 for the Landau-Ginzburg parameter. As is 
known from conventional superconductors, if ft > 1/V2, 
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distant vortices repel each other leading to type-II be- 
havior. This is the standard picture of the proton super- 
conductor in neutron stars that is widely accepted in the 
astrophysics community. 

However, the standard procedure described above is 
inherently flawed since the system exhibits an approxi- 
mate U(2) symmetry, and therefore the couplings are 
approximately equal an m a 22 ~ 0-12- This makes the 
mixing matrix M nearly degenerate. The general solu- 
tion to Eq. (O is: 

(g) = ^ CiK {^ir) v ; (22) 

where v% and Vi are the eigenvalues and eigenvectors 
of matrix M, and Cj are constants to be calculated by 
matching to the solution of the original nonlinear equa- 
tions of motion. We would like to introduce the parame- 
ter e (which measures the asymmetry between the proton 
and neutron Cooper pairs) defined in the following way: 

Sd 

e = (011022 - alz)/ a ii — 2 — ■ ( 23 ) 

a 

We should remark here that our qualitative results do 
not depend on the value of 7 = ni/n 2 . Indeed, no matter 
what ni and n 2 are, the mixing matrix is still singular 
in the limit e — > 0. Hence, we still get one eigenvalue 
which vanishes when e — * 0. So, the only crucial as- 
sumption is e <C 1. However, to simplify our formula 
for the eigenvalues in what follows we assume a specific 
value of 7 = n\jn 2 C 1. It simply allows our results 
to be expressed in a more transparent way. In the limit 
7 = ri\Jn% <C 1 and e = 25a ja -C 1 one can estimate the 
eigenvalues and eigenvectors of the matrix M as: 

-1 - p, v. - (^) , (24) 

,2 - ^ v 2 -(;). (25) 

The physical meaning of solution (1221) is simple: there are 
two modes in our two component system. The first mode 
describes fluctuations of relative density (concentration) 
of two components and the second mode describes fluctu- 
ations of overall density of two components. Notice that 
v\ <C V2, and hence the overall density mode has a much 
smaller correlation length than the concentration mode. 
Therefore, far from the vortex core, the contribution of 
the overall density mode can be neglected, and one can 
write: 

(F) (v ► 00) — CiKoiVTer/O ■ ('J) (26) 

The most important result of the above discussion is that 
the distance scale over which the proton and neutron con- 
densates tend to their vacuum expectation values near a 
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FIG. 1: In this figure we show the functions /(f), G(r)/j, and 
a(f) (defined in Eqs. 1131161 ') as a function of the dimension- 
less radial coordinate r — r/£. The dotted line corresponds 
to a(f), the solid line approaching 1 at large f corresponds to 
/(f), and the solid line approaching at large f corresponds 
to G(f)/7. 



proton vortex is of order - the correlation length of 
the concentration mode. Since e -C 1, this distance scale 
can be much larger than the proton correlation length £, 
which is typically assumed to be the radius of the proton 
vortex core. The appearance of the concentration mode 
is not surprising since we presented some arguments (be- 
fore calculations) supporting the picture that the neutron 
condensate will increase its magnitude slightly in the vor- 
tex core, while the proton condensate will decrease its 
magnitude to in the core center. We note that in the 
limit e — > the size of the proton vortex core becomes 
infinitely large, and the vortex is thereby destroyed. This 
is in accordance with the topological arguments, which 
state that if the U(2) symmetry were exact with e = 0, 
and it is spontaneously broken to U(l), there will be no 
stable vortices in the system. 

We have also verified the above results numerically by 
solving the equations of motions 11UI11I12|I with a par- 
ticular choice of the approximately U(2) symmetric in- 
teraction potential V. Our numerical results support the 
analytical calculations given above. Namely, we find that 
the magnitude of the neutron condensate is slightly in- 
creased in the vortex core, the radius of the magnetic flux 
tube is of order A and the radius of the proton vortex core 
is of order £/t/e. In Fig. \I\we show the numerical solu- 
tion of the profiles of the proton vortex (/(f)), neutron 
condensate (G(f)/j), and a(r) (related to the gauge field 
through Eq. Ultty l as a function of the dimensionless ra- 
dial coordinate f = r/£, where £ is the coherence length 
(nj. We have used « = 3, n 1 /n 2 = 0.05, and e = 0.02 in 
this numerical solution. 



III. VORTEX- VORTEX INTERACTION 

Now that we know the approximate solution for the 
proton vortex, we will proceed to look at the interac- 
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tion between two widely separated proton vortices. If 
the interaction between two vortices is repulsive, it is en- 
ergetically favorable for the superconductor to organize 
an Abrikosov vortex lattice with each vortex carrying 
a single magnetic flux quantum. As the magnetic field 
is increased, more vortices will appear in the material. 
This is classic type-II behavior. If the interaction be- 
tween two vortices is attractive, it is energetically favor- 
able for n vortices to coalesce and form a vortex of wind- 
ing number n. This is type-I behavior. Typically, the 
Landau-Ginzburg parameter k = A/£ is introduced. In 
a conventional superconductor, if k < 1/ 1 \[2 then the su- 
perconductor is type-I and vortices attract. If k > 1/V2 
then vortices repel each other and the superconductor is 
type-II. As mentioned above, the typical value for a neu- 
tron star is n ~ 3, so one could naively expect that the 
proton superfluid is a type-II superconductor. 

The case we considered in the previous section has one 
new element, e which was not present in the standard 
type-I/II classification. However, we expect that analo- 
gous classification should remain in effect for the proton 
vortices described above. In such an analysis the coher- 
ence length £ should be replaced by the actual size of the 
proton vortices 5 ~ C/v 7 ^- Therefore, we will define a 
new Landau-Ginzburg parameter for our case, 



_ A _ 



-A 



(27) 



We expect type-I behavior with vortices attracting each 
other if K np <C 1 and type-II behavior if n np 1. For 
relatively small e such an argument would immediately 
suggest that for the typical parameters of the neutron 
stars type-I superconductivity is realized (rather than the 
naively assumed type-II superconductivity) . In what fol- 
lows we present several different arguments supporting 
this claim. 



Calculating in the intervortex potential: 
method I 



eliminating the phase of ipi field in favor of the longitu- 
dinal component of A, to produce a non-interacting free 
energy T free . 



T 



free 



2m, 



■(m(VF) 2 + n 2 (VG!) 2 ) 



l (( VxAf + lA 2 ) + 2«^ 



+ 4(a - 8a)nin 2 FG + 2anjG 2 



(28) 



Following 0, 0| , we must also add the source terms for 
each field to model the vortices: 



J d 2 r(pF + tG + j • A) 



(29) 



The solutions to the equations of motion following from 
Tfree coupled to the sources can be obtained in the same 
manner as was done in 14). The equations of motion 
resulting from Tf Tee + Source are: 



(V 2 - M) 



~h? 



p/m 
r/n 2 



V 2 A 



A 2 



A = 47TJ 



(30) 
(31) 



where M is the same mixing matrix given in Eq. I|19|l . 
Since we are interested in the asymptotic behavior, we 
will choose the the first eigenvalue/eigenvector solution 
(|24|l that diagonalized M. Therefore, the asymptotic 
field solutions are given by Eqs. (|2UI26I) : 



(32) 



F ~ -G/ 7 ~ -C 1 ^o(V2er/£) 

he 



A ~ —C A Kx(r/\) e. 
q\ 

Following [l^, we require that (|32|l are solution of the 
equations of motion l|30|l and l|31|l . The source solutions 
can immediately be written down when we compare Eqs. 
(|30|l and (|31|l with the following Bessel equations: 



In order to make this qualitative discussion more con- 
crete, we will present three different calculations sup- 
porting our claim that for the typical parameters of a 
neutron star the proton superconductor may be type-I 
rather than type-II. First of all, we follow the method 
suggested originally in 0] to calculate the force between 
two widely separated vortices. The methods of 0] were 
subsequently applied in 0| to the case similar to ours, 
the interaction of two widely separated vortices that have 
nontrivial core structure. In these papers, the force be- 
tween two widely separated vortices is calculated by using 
a linearized theory with point sources added at the loca- 
tion of the vortices. The point sources in the linearized 
theory are chosen to produce fields matching the long 
distance asymptotics of the original theory. Therefore, 
we expand the free energy up to quadratic order in 
the fields F, G, and A introduced in the previous section 



(V 2 - p 2 )K (fix) = -2tt<5(x) 



(33) 



(V 2 - ^K^x) 

x fj, 



eV(x) (34) 



The source solutions that solve the equations of motion 
along with (I32|) are: 



5® 

he 

- C A Vx(S 2 (r)z) 
2q 



(35) 



Since the equations of motion corresponding to J- free 
coupled to T 'source are linear, the two- vortex solution can 
be considered as the sum of two single vortices at posi- 
tions i - ! and r 2 . To calculate the vortex- vortex interac- 
tion energy, we use ansatz (F, G, A) = (Fi + F 2 , G\ + 
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G 2 ,Ai + A 2 ) (p,r,j) = (px + p 2 ,n + t 2 Ji + j 2 ) in 
Ffree + ^source an d subtract off the energy of each iso- 
lated vortex. The notation 1,2 indicates that these are 
functions of r — 1*1,2 ■ Using the equations of motion 
(|30I31|I the interaction energy can be written as 



Til 



<rr(ji • A 2 + P1F2 + nG 2 ) 



(36) 



Substituting the asymptotic field solutions (|32I35|) into 
the interaction energy given above, the integration can be 
done as in 01 to obtain the following expression for the 
interaction energy per unit vortex length of two widely 
separated parallel vortices: 



U(d) ~ ^±(C%K (d/X) 
m c 

- C*(l + O( 7 ))K (V2ld/t)) 



(37) 



where d — \v\ — r 2 | — > 00 is the separation between the 
two vortices. We see that if the first term in U dominates 
as d — > 00 then the potential is repulsive, otherwise, if 
the second term dominates the potential is attractive. 
In other words, if < 1/V%, then vortices attract 

each other and the superconductor is type-I; otherwise, 
vortices repel each other and the superconductor is type- 
II. This confirms our original qualitative argument that 
K 7i P = A/<5 = y/e\/£ should be considered as an effec- 
tive Landau- Ginzburg parameter, which determines the 
boundary between the type-I and type-II proton super- 
conductivity. In terms of the parameters of our theory, 
we have 



A 
~5 



X 



6 a 



(38) 



In this case we see that the type-I/II behavior is con- 
trolled in part by the degree of symmetry breaking (pro- 
portional to Sa). Our numerical estimates (see conclu- 
sion) suggest that n np < l/\/2, and therefore, the system 
is a type-I superconductor. 



B. Calculating the intervortex potential: method II 

Due to the importance and far reaching consequences 
of this result, we have also calculated the vortex- vortex 
interaction energy in a more direct way following |l5j| . 
where the introduction of the auxiliary sources is com- 
pletely avoided. This method, very different in nature, 
has produced the same result (|37|l as the above proce- 
dure, therefore confirming our picture. 



We proceed as follows. First, we again eliminate the 
phase of tp\ with an appropriate gauge transformation. 
Now, recalling our definitions l|13l 1161) , let (F, G, A) = 
(Fi , G\ , Ai ) be the exact fields of a single vortex located 
at n. Also let (F, G, A) = (F a + F 2 , G x + G 2 , Ai + A 2 ) 
be the exact fields produced by two vortices at locations 
ri and r 2 . Note that subscript 2 does not generally refer 
here to functions of r— r 2 as it did in the previous subsec- 
tion: here (F 2 ,G 2 , A 2 ) are just corrections to the fields 
of the single vortex at ri . When r is far from the core of 
vortex 1, Fi, Gi, Ai are small and the equations of mo- 
tion can be linearized to yield familiar asymptotics l|32f) . 
Moreover, when r is far from the cores of both vortices, 
Fi, Gi,A\ are small for both i = 1, 2 and we obtain: 



-Gih -dtfotV&lr-nl/O, 
he 

A ; ~ -^Mxdr-ril/A) 6. 
qX 



(39) 



To calculate the interaction energy of the two vortices, 
let us divide the space into two cells Ti and T 2 , which 
contain the centers of vortices 1 and 2 respectively (if we 
had more than two vortices in our problem these would 
be the Wigner-Seitz cells of the vortex lattice). For sim- 
plicity, we will take the boundary between two cells to be 
the line perpendicular to and passing through the mid- 
dle of the line-segment joining the cores of two vortices. 
This boundary then acts as an axis of symmetry for the 
problem: the total energy of two interacting vortices is 
just twice the energy of each cell. The vortex-vortex in- 
teraction energy is then: 

T int = T[F 1 +F 2 ,G 1 +G 2 ,A 1 +A 2 }-2T[F 1 ,G 1 ,A 1 } 
= 2 / d 2 r(£[F 1 +F 2 ,G 1 + G 2 ,A 1 +A 2 ] 

- £[F 1 ,G 1 ,A 1 ])-2 f d 2 r£[F 1 ,G 1 ,A 1 ] (40) 

where £ — dT jdr 2 is the free energy per unit volume and 
subscripts Ti indicate integration over the corresponding 
cells. We assume the vortex separation to be large, there- 
fore the fields F\,Gi,A\ inside cell T 2 will be small, so 
we can expand the integral over the cell T 2 in 14011 to 
second order in Fi,Gi,Ai. Similarly, the corrections to 
fields Fi,Gi,Ai inside cell T\ will be small so we can 
expand the integral over the cell Ti in (|40() to second 
order in corrections F 2 ,G 2 , A 2 . The resulting expression 
can be simplified by using equations of motion l|10f) , pif) > 
{El to obtain: 
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J~ it 



JTi 2m e 

dS ^ 2m c 



F 2 {2VF 1 + VF 2 



2m r 



G 2 (2VGi + VG 2 



<Stt 



A 2 x (2V x Ai + V x A 2 )) 



, G1VG1 + — A x x (V x AO) 



(41) 



Here the integrals are over the boundary of cells T\ and 
T 2 . Since this boundary is far away from either vor- 
tex center, we can use the asymptotic expressions (|39f) 
for the fields (Fj,Gi,Ai) to explicitly calculate the sur- 
face integrals in 1)41 [I. We note that by symmetry of the 
asymptotic solution on the boundary, the second integral 
in (|41|l cancels with the part of the first integral to yield: 



J' it 



\l dS • ( 



^1 Wl + ^G 2 V Gl 



— A 2 x (V x AO) 

47T 



(42) 



Substituting asymptotic solutions into the above, we 
find the vortex-vortex interaction energy per unit length 
to be: 



U(d) 



(CiK (d/X) 



by Eq. QJ) over a superconducting sample. For the su- 



perconducting state where (\ipi 



ni, (\ip2 



ri2, 



and B — (Meissner effect), the Gibbs free energy is 



9s{H,T) 



.IL 

2a 



5a 



\2a) 



(45) 



where we expressed the result in terms of the parameters 
of Eq. 0. We also replaced the densities (m + n 2 ) — > 
/x/a and (n 2 — ni) — > 25fi/5a in terms of the same param- 
eters by neglecting small factors ~ da/ a. For the normal 
state, we have (iV'il 2 ) = 0, (\ip2 1 2 ) = (m + <5pt) /a, and 
B = H. The Gibbs free energy is: 

^! 

' 8tt 



9n(H,T) 



2a 



(46) 



As H c is defined as the point at which g s (H c ) = g n (H c ), 
we can solve for the critical field. Equating the free en- 
ergies of the normal and superconducting state, we find 



- C 2 (l + 0( 7 ))ifo(V27d/0) (43) 
in accordance with our previous result 



IV. CRITICAL MAGNETIC FIELDS 

Our third check of the main result that for relatively 
small e the superconductor in the neutron stars may be, 
in fact, type-I is based on the calculation of the criti- 
cal magnetic fields. Usually one calculates the critical 
magnetic fields H c and i? c2 . These are the physically 
meaningful fields above which the superconductivity is 
destroyed in type-I and type-II superconductors respec- 
tively. If H c > H C 2 then the superconductor is type-I, 
otherwise, the superconductor is type-II. 

First, we will calculate the critical magnetic field H c . 
This is defined as the point at which the Gibbs free energy 
of the normal phase is equal to the Gibbs free energy 
of the superconducting phase. In other words, as the 
external magnetic field H is increased above H c , it is 
energetically favorable for the superconducting state to 
be destroyed macroscopically. The Gibbs free energy in 
the presence of an external magnetic field H is: 



g(H,T) = f(B,T)-™ 



(44) 



where H is the external magnetic field, B is the magnetic 
induction, and T denotes temperature. The quantity 
f(B,T) is the integrand of the free energy density (given 



11. =N — (ii-^f 
2a da 



niV87T(5a, (47) 



where at the final stage we used the equation for n\ in the 
superconducting phase expressed in terms of the original 
parameters J3J|. 

Now we will proceed to calculate i/ c2 . This is the crit- 
ical magnetic field below which it becomes energetically 
favorable for a microscopic region of the superconduct- 
ing state to be nucleated, with the normal state occurring 
everywhere else in space. In order to calculate i? c2 , we 
follow the standard procedure and linearize the equations 
of motion for tpi about the normal state with (|'0i| 2 ) = 
and (|^ 2 | 2 ) = (p + 5/j,)/a. The linearized equation of 
motion reads, 



2m r 



q \* 
-iv - — A J Vi = wipi, 

uj = (fj, + Sfj,) 2<5/i. 



(48) 
(49) 



This is simply a Schrodinger equation for a particle in 
a magnetic field, with an energy of u). This is a stan- 
dard quantum mechanics problem and we can immedi- 
ately write down the solution. The first Landau level is 
the ground state energy of cq(H) = h\q\H/2m c c. There- 
fore, if u) < 6q, then only the trivial solution with ip\ = 
is possible. The critical field H C 2 is defined as the point 
at which u) = eo(H C 2). This is given as 

2m c c,, . .5a nc . Am c c . 

H c2 = —^-[(fi + Sfi) 2fy-— -jfom 50 

n\q\ a n\q\ 
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Now that we have the critical fields H c and H C 2 in hand, 
we can compare the two in order to determine the typc- 
I/II nature. If H c < H C 2 this means that it is energet- 
ically favorable for microscopic regions of the supercon- 
ducting state to be nucleated as H is decreased. This 
is type-II behavior, and this nucleation manifests itself 
in the form of an vortex lattice. If H c > H C 2, then it 
is energetically favorable for macroscopic regions of the 
superconducting state to be present as H is decreased. 
This is a type-I superconductor and the superconducting 
state persists everywhere in the material when H < H c , 
as opposed to a type-II superconductor where it is local- 
ized in regions of space in between the vortices. From 
Eqs. P47I49I50|) we see that 



H c2 ^ ^ TOcC 



1 6a 



y/nh q 



np • 



(51) 



This agrees with the parametrical behavior given in Eq. 
(|38|l obtained from the vortex interaction calculation of 
the previous section. To estimate H c numerically, it is 
convenient to represent H c as 



H r . 



2ir\£ 



</>o 



2-kKc 

q 



2 x 10 7 G • cm 2 , (52) 



where ipo is the quantum of the fundamental flux. If we 
substitute A = 80 fm and £ = 30 fm (typical values) in 
the expression for the critical magnetic field (|47|l . H c is 
estimated to be the H c ~ 10 14 G, which is smaller than 
the "naive" estimate by a factor of y/Sa/a ~ 10 _1 . It is 
quite amazing that very different calculations of the crit- 
ical magnetic fields (|51|1 lead exactly to the same conclu- 
sion which was derived from analysis of the vortex- vortex 
interaction (|38[) . 



V. CONCLUSION 

In this paper we have demonstrated using various cal- 
culations that the proton superfluid present inside a neu- 
tron star may in fact be a type-I superconductor Q . This 
supports the observation made my Link Q that the con- 
ventional picture of type-II superconductivity may be in- 
consistent with the observations of long period precession 
in isolated pulsars @. The most important consequence 
of this paper is that whether the proton superconductor 
is type-I or type-II depends strongly on the magnitude of 
the SU(2) asymmetry parameter e. Specifically, we find 
that the superconductor is type-I when K np = -y/eA/f; < 
l/y/2, and type-II otherwise. This result is quite generic, 
and not very sensitive to the specific details of the inter- 
action potential V. In particular, when e — > the super- 
conductor is type-I. The parameter e is not known pre- 
cisely; the corresponding microscopical calculation would 
require the analysis of the scattering lengths of Cooper 
pairs for different species. We can roughly estimate this 
parameter as being related to the original SU(2) isospin 
symmetry breaking e ~ (m n — m p )/ra n ~ 10~ 2 . If we 



assume a typical value for A/£ ~ 3 and e ~ 10~ 2 , we 
estimate K np — ~ 0-3 < l/v2> which corresponds 

to a type-I superconductor. From these crude estimates, 
we see that it is very likely that neutron stars are type- 
I superconductors with the superconducting region de- 
void of any magnetic flux, as was originally suggested in 
to resolve the inconsistency with observations of long 
period precession j3] in isolated pulsars. Many mecha- 
nisms have been proposed to explain glitches [lTj ■ If the 
proton superfluid exhibits type-I superconductivity then 
some explanations of glitches [l8j that assume type-II 
superconductivity would have to be reconsidered, as sug- 
gested in Ref. pj . It might be interesting to consider how 
the presence of a nonzero proton condensate affects the 
characteristics of the neutron vortices that carry angu- 
lar momentum. In particular, as we have demonstrated 
for magnetic flux tubes in this paper, the neutron vor- 
tices might have an enhanced proton superfluid density 
inside, as well as a coherence length (the approximate 
size of the vortex) which is much larger than originally 
expected. In this case, the pinning force that is related 
to the size of the vortex core |l9j could be very different 
from the simple estimates when the strong interaction 
between the neutron and proton Cooper pairs is ignored. 

If the core is indeed a type-I superconductor, the mag- 
netic field could exist in macroscopically large regions 
where there are alternating domains of superconducting 
(type-I) matter and normal matter. In this case, neutron 
stars could have long period precession. Such a structure 
follows from few different arguments. First of all, as has 
been estimated 20], it takes a very long time to expel a 
typical magnetic flux from the neutron star core. There- 
fore, if the magnetic field existed before the neutron star 
became a type-I superconductor, it is likely that mag- 
netic field will remain there. Another argument suggest- 
ing the same outcome follows from the fact that topol- 
ogy (magnetic helicity) is frozen in the environment with 
high conductivity; therefore, the magnetic field must re- 
main in the bulk of the neutron star. The last argument 
supporting the same picture is due to Landau |2l| who 
argued that if a body of arbitrary shape (being a type-I 
superconductor) is placed under influence of the external 
magnetic field with strength H < H c , then the magnetic 
field in some parts of the body may reach the critical 
value H c , while in other parts of the body it may remain 
smaller than the critical value, H < H c . In this case, 
a domain structure will be formed, similar to ferromag- 
netic systems. Specifically, on a macroscopic distance 
scale, the magnetic flux must be embedded in the super- 
conductor. This would mean that the superconductor is 
in an intermediate state as opposed to the vortex state 
of the type-II superconductor, which was assumed to be 
realized up to this point. The intermediate state is char- 
acterized by alternating domains of superconducting and 
normal matter. The superconducting domains will then 
exhibit the Meissner effect, while the normal domains will 
carry the required magnetic flux. The pattern of these 
domains is usually strongly related to the geometry of 
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the problem. The simplest geometry, originally consid- 
ered by Landau |2l| , is a laminar structure of alternating 
superconducting and normal layers. However, it has been 
also argued that due to the geometry of the neutron 
vortex lattice, the normal proton domains will be in the 
shape of cylinders concentric with the rotational neutron 
vortices. 

While some precise calculations are required for under- 
standing of the magnetic structure in this case, one can 
give some simple estimation of the size of the domains 
using the calculations Landau presented for a different 
geometry. His formula [2l| suggests that the typical size 
of a domain is 

a ~ 10V#A, (53) 

where R is a typical external size identified with a neu- 
tron star core (R ~ 10 km), while A is a typical width of 
the domain wall separating normal and superconducting 
states. We estimate A ~ 5 = £/y/e as the largest micro- 
scopical scale of the problem. Numerically, a ~ 10 _1 cm 
which implies that a typical domain with size of order 
~ 10 _1 cm can accommodate ~ 10 neutron vortices sep- 
arated by a distance ~ 10 -2 cm. 

The consequences of this picture still remain to be ex- 
plored. In particular, it might be of importance in the 
explanation of glitches. It may be also important in anal- 
ysis of the cooling properties of the neutron stars. 



It would be very interesting to test the ideas outlined 
in this paper by doing laboratory experiment in the spirit 
of the Cosmology in the Laboratory (COSLAB) program. 
In particular, it would be interesting to find a condensed 
matter system (high T c superconductor?) where the very 
interesting feature discussed in this paper can be tested. 
Namely, the core of the vortices and their interactions 
could be very different from the simple estimates. This 
might happen if there is a condensate of another field 
with energy scales almost degenerate with energies de- 
termined by the Landau- Ginsburg complex ip field de- 
scribing a superconductor. Over the last few years several 
experiments have been done to test ideas drawn from cos- 
mology and astrophysics (see Ref. Q and the web page 
of the latest COSLAB meeting for further details). 
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